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PERFORMANCE OF TWIN GYRO ATTITUDE CONTROL SYSTEM 
INCLUDING PASSIVE COMPENSATION AND 
NONLINEAR CONTROL LAW 

INTRODUCTION 


Newton's second law, the Principle of Conservation of Angular Momentum, 
forms the foundation upon which all attitude control systems are based. These 
systems are categorized as being either active, passive, or semipassive. 

An active control system is defined as one that requires some form of 
onboard power. These are generally closed-loop systems requiring some form 
of sensor to provide the feedback signal [ 1] ; that is, instrument gyro, sun sen- 
sor, horizon sensor, or star tracker. These types of systems are employed for 
missions requiring high pointing accuracy or for performing complex maneuvers. 
The torque -producing devices for such systems are reaction control jets (RCJ's) , 
reaction wheels, or control moment gyros (CMG's) . 

Passive control systems are defined as those systems that do not require 
onboard power, but rather they make use of such natural phenomena as solar 
radiation, magnetic gradients, and gravity gradient [2-4]. These are open-loop 
systems with very large time constants. Systems of this type have been used on 
weather satellites where long mission life and high reliability are key factors. 

In some cases it is possible to replace active components with passive 
components, thereby forming a hybrid system referred to as "semipassive." 
These systems take advantage of the power savings and high reliability of passive 
components and yet retain the advantages of the active systems. 

For vehicles requiring very accurate and continuous control, the CMG 
is the superior device from the standpoints of power, accuracy, response, and 
simplicity [5]. 

A CMG consists basically of a rotor with a certain constant magnitude of 
angular momentum and a mechanism used to precess the rotor. The configuration 
of this precession mechanism determines the cross coupling characteristics of 
the generated control torques. There are three basic configurations of preces- 
sion mechanisms referred to as one degree of freedom, two degrees of freedom, 
and the twin gyro (Fig. i) . 



The one-degree -of -freedom CMG, Figure la, consists of a rotor mounted 
on a gimbaled shaft driven by a torque motor. The developed control torque lies 
in a plane perpendicular to the gimbal shaft, thus producing two-axis cross cou- 
pling. This control torque is transmitted through the gimbal axis bearings to the 
vehicle. A control system composed of three one-degree-of-freedom CMG's was 
discussed and analyzed by White and Hansen [ 6] . These systems were shown to 
be severely cross coupled and required a very complex control computer to 
decouple the control torques so that three-axis control could be attained. A 
failure of any one CMG would cause the loss of three-axis control. 

The two-degrees-of-freedom CMG consists of a rotor mounted so that it 
is capable of being precessed about two axes. This is accomplished by the double- 
gimbal precession mechanism shown in Figure lb. This mechanism requires 
two torque motors, one for each of the gimbal axes, which must be capable of 
transmitting the control torques to the vehicle. This type of precession mecha- 
nism produces three-axis cross coupling of the control torque. A minimum of two 
two-degrees-of-freedom CMG's are required to obtain three -axis attitude con- 
trol. This system is also severely cross coupled and requires a very complex 
control computer for its implementation [7, 8] . Unlike the one-degree-of- 
freedom CMG system, this system is 33 percent redundant with respect to its 
torque motors; that is, should one motor fail, this system would still be able to 
provide three -axis control. 

The third configuration, the twin gyro shown in Figure lc, consists of two 
one-degree-of-freedom CMG's geared together so that their angular momentum 
vectors always remain in a position that is the mirror image of the other. This 
configuration, like the one degree of freedom, requires only one torque motor 
and permits the control torques to be transmitted through Ihe gimbal shaft 
bearings to the vehicle. The major advantage of this configuration over the 
previous two is that it eliminates the first-order cross coupling terms, thereby 
forming a single -axis control device. A minimum of three twin CMG's are 
required for three-axis attitude control [9] . This control system possesses 
only second-order cross coupling effects and therefore requires the simplest 
control computer. While it is true that the weight and power requirements are 
greater for this system than for either of the two systems mentioned previously, 
this system does possess 100 percent redundancy with respect to its rotors; that 
is, three-axis control can be maintained with one rotor failure for each twin 
gyro. However, with the loss of a rotor the twin gyro is transformed into a 
single-degree-of -freedom gyro; thus, an alternate control law must be imple- 
mented by the control computer to facilitate the two-axis cross coupling present. 

It is the opinion of H.B. Kennedy [ 10] that the reliability of the twin CMG is 
greater than that of the two -degree -of-freedom CMG because of the reduced 
number of gimbals and gyro bearings . 


2 





b. Two Degrees of Freedom 



c . Twin Gyro 


Figure 1. Control moment gyros. 
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Throughout this paper the complete assemblage of the twin CMG will be 
referred to as being a gyro platform. When these platforms are fastened rigidly 
to the vehicle structure, they are referred to as being rigid. When the platforms 
are fastened by means of a spring-damper mechanism, they are referred to as 
being compensated. 

DERIVATION OF CONTROL EQUATIONS FOR A VEHICLE 
WITH RIGID PLATFORMS 


Consider a three-axis attitude control system which is composed of 
three twin CMG platforms, as shown in Figure 2. The assumptions made for 
the derivation of the equations of motion for this system are as follows: 

1. The gyro platforms are rigidly fastened to the vehicle structure; 
that is, the vehicle-fixed coordinate system corresponds to the gyro-platform- 
fixed coordinate system. 

2. The mass moments of inertia of each rotor and each of their gimbal 
structures are equal. 

3. The magnitudes of the rotor's angular momentum vectors are equal 
and constant. 

4. No mass unbalance exists about the gimbal axis; that is, the center of 
mass of the gyros lies on the gimbal axis. 

The three CMG platforms shown in Figure 2 are designated as X, Y, and 
Z, corresponding to the vehicle's X-, Y-, and Z -axes, respectively. These 
platforms are aligned so that the developed control torques act along the vehicle's 
principal axis. 

A vector schematic of the X-axis platform is shown in Figure 3 . The 

angular momentum vector of each gyro is designated H and H . The 

X 1 x 2 

first subscript associates the angular momentum vector with a particular plat- 
form; the second subscript indicates the gimbal axis, axis i being the torquer- 
driven shaft. The angular momentum vectors, written in terms of the vehicle- 
fixed coordinates, are 
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( 1 ) 


H = H cos 6 i + H sin 5 £ 
Xj v x t v 


— A A 

H = - H cos 6 i + H sin 6k. (2) 

x 2 x 2 J v x 2 v 


From classical dynamics the torque developed by precession of an angular mo- 
mentum vector, with respect to an inertial reference, is given as [ ii] 


r = 


dH 

dt 


Ref 


+ co x H 


(3) 


Note that equation (3) yields the torques acting on the gyro rotors. To 
obtain the torques acting on the vehicle, one must take the negative of equation 
(3) , since the vehicle control torques are actually the reaction torques of those 
described by equation (3) . Therefore 


T 

Control 


dH 

dt 


Ref 


- CO X H 


(4) 


The first term in equation (4), 


dH 

dt 


Ref’ 


gives the torques produced by the gyro's 


precession mechanism, while the second term, co x H, is produced by the angu- 
lar velocity (co) of the vehicle. The angular velocity of the vehicle is defined 
as 


— ► A 

co = r i 


A 

Pj 


A 

+ qk 
v v 


(5) 


where r, p, and q are the angular velocity components in the vehicle -fixed 
coordinate system. The torques acting on the rotors are then found by substi- 
tuting equations ( i) , ( 2) , and ( 5) into equation ( 3) , yielding 

T = H (psin 6 - qcos 6) i - H (r + 6) sin 6 j +H (r + 6) cos 
Xj Xj n V X! J V Xi v 

( 6 ) 


r = H (psin 6 + qcos 6) 1 - H ( r - 6) sin 6 j - H (r - 6) cos 6k 
x 2 x 2 K v x 2 J v x 2 v 

(7) 
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A free body sketch of the X-axis platform is shown in Figure 4. The torque 
components of equations (6) and (7) are shown acting on their respective rotors. 
The inverse of these rotor torques is then shown acting on the gimbal shaft's 
structure. The result of these rotor torques, as reflected on the number one 
gimbal shaft, is obtained by adding their respective components. Note that the 
torques about the Y- and Z-axis add directly but, because of the gear configura- 
tion, a positive torque about the X-axis of the number two gimbal shaft will be 
reflected onto the number one shaft as a negative; thus 


T V =T V -T V =-H (psin 6 - qcos 6)+H (psin 6 + qcos 6) 

XX XjX XgX x t Xg 


T V = T V + T V = H (r + 6) sin 6 + H ( r - 6) sin 6 
xy Xtf X 2 y X! x 2 


T V = T V +T V = -H (r + 6) cos 6 + H ( r - 6) cos 6 

XZ XjZ X 2 Z Xj x 2 


( 8 ) 

(9) 

( 10 ) 


The second subscript indicates the component of the torque with respect to the 
first subscript's platform. The superscript v indicates the vehicle coordinate 
system. 

Since it was assumed that the magnitude of the rotor's angular momentum 
vectors were equal, 


|H 


xi' 


= I H 


x 2 


= IH 


Equations (8), (9), and (10) are reduced to 


T V = 2 H q cos 6 

XX X 


(11) 

T V = 2H r sin 6 
xy x 


(12) 

T V = - 2H 6 cos 6 

XZ X 

• 

(13) 
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V V 

The torques T and T are transmitted through the gimbal axis bearings to 
xy y ^ 

the vehicle, while T acts upon the gimbal shaft. 

XX 

At this time it is necessary to derive the equation of motion for the gim- 
bal shaft. This is accomplished by summing torques about the X-axis of the 
gimbal shaft and equating them to the Euler equation of motion. A generalized 
form of Euler's equations of motion for a rigid body whose XYZ axes are aligned 
with the principal axes and whose origin is coincident with the center of mass is 
given as [ 11] 

ST x =I x t + pq(I z“ I y ) U4) 

ST = I p + rq(I - I ) (15) 

y T X z 

ST = I q + rp(I - I ) (16) 

z z y x 

Summing the torques shown in Figure 4 yields 

T - D 6 + 2H q cos 6 = I g 6 . (17) 

mx x x xx 


T^^ is the torque applied by the gimbal shaft torque motor, and is a 

rate damping torque created by an eddy current damper, whose damping coef- 
ficient D has the units • The moments of inertia of the gimbal shaft 

g g g 

associated with the X-axis platform are denoted as I , I , and I . The 

xx xy xz 

first subscript defines the platform, and the second subscript defines the direc- 
tion in that platform. The torque described by the term pq^I^ - has been 

neglected in equation ( 17) as being insignificant in comparison with the other 
terms. This term's insignificance is a result of both the low vehicle rates and 
the very small mass moments of inertia of the gimbal shaft. 
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The torque acting on the vehicle in the X-direction can now be seen to be 
the reaction torques of the torque motor and the eddy current damper. Hence the 
control torques applied to the vehicle are 


' v = -T + D 6 

(18) 

xx mx x 


’ V = 2H r sin 6 

( 19 ) 

xy x 


V * 

' = -2H 6 cos 6 

(20) 

xz X 



These torques are shown superimposed on the X-axis platform ( Fig. 5) . 

The control torques developed by the Y- and Z-axis platforms are 
derived in the same manner. Vector diagrams of these platforms are shown in 
Figures 6 and 7 where jj. and y are the respective gimbal angles. The equa- 
tions of motion for their respective gimbal shafts are 


T -D/I+2H rcosn = I g |I (21) 

my y y yy 


T - D v + 2H p cos y = I ® y 
mz z z zz 


( 22 ) 


The developed control torques are as follows: 


Y-axis platform 


, v ... • 

1 = -2H n cos M 

(23) 

yx y 


' v = -T + D A 

(24) 

yy my y 


’ = 2H p sin ju , 

(25) 

yz y 
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( 26 ) 


T V = 2H q sin y 
zx z 


T V = -2H y cos y ( 2 ?) 

zy z 


T V = -T + D y . (28) 

zz mz z 


These torques are also shown on their respective platforms ( Fig. 5) . 


Z VEHICLE 



Figure 5. Schematic of the control torques generated by each gyro platform. 
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Z INERTIAL 
REFERENCE 



Figure 6. Vectorial illustration of Y-axis platform. 

The equations of motion for the vehicle are obtained by summing the 
torques shown on Figure 5 and equating them to Euler's equations of motion. 
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External disturbance torques (T ) are assumed to oppose positive gyro con- 
tol torques. 
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X-axis 


2 H z q sin y - 2H L cos \x - 1^+ D* + - fi + pq (ij - ij) ( 29) 

X 

Y-axis 

2 H x r sin « - 2H z r cos T - T my + D/ + T - £ P + *<l(£ ' '!) ( 30) 

Z-axis 

2H p sin m - 2H 6 cos 6-T + D „ T+ T , =I V q+ rpfl V - I V ) ( 31) 

y x mz z ext z V y x/ 

z 

The terms I V , I V , and I V represent the mass moments of inertia of the 
x y z 

vehicle about its X-, Y-, and Z-axes, respectively. Figure 8 is a block diagram 

of the three-axis control system. Because the torque motor's time constant is 

so very small, a constant, K , will be used for its transfer function. The 

m 

vehicle rates r, p, and q are transformed into Euler rates using the following 
transformation equations [11]: 


^ = (p sin <f> + q cos 4>) sec 0 

9 = p cos 4> - q sin 4> 

• • 

<t> = r + ip sin 9 


(32) 

(33) 

(34) 


The integral of these Euler rates is used as the feedback to close the loop about 
the system. A constant gain is applied to each of the feedback signals. 

From the generalized vehicle equations of motion [equations (29) , (30) , 
and (31) ] , it can be seen that if the angular velocity of the vehicle is initially 
zero and an attitude maneuver is commanded about any one of the vehicle's prin- 
cipal axes, the equation of motion about that axis can be closely approximated as 

X-axis 

T - 2H (i cos M = I V r ( 35) 

ext y x 

x 

Y-axis 

T - 2H 7 cos y = I V p (36) 

ext z y 
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Figure 8. Three-axis rigid platform control system. 
















Z-axis 


T - 2H 6 cos $ = I V q 
ext x z 

z 


( 37 ) 


These equations are defined as the uncoupled, single-axis equations of motion 
for the vehicle. 


The simplifications made are justified when one considers that the mag- 
nitude of the reaction torque of the torque motor and eddy current damper are 
extremely small (two orders of magnitude less than the principal control torque 
for the system to be discussed later) , and therefore its effect on the large mo- 
ment of inertia of the vehicle is negligible. Figure 9 is a block diagram of the 
uncoupled Z-axis control system. 


This is a dual input system; one is a vehicle attitude command, 0 , and 

c 

the other is an external disturbance torque. Assuming that 


T ex t (s) ' 0 • 

the forward- and open-loop transfer functions are, respectively, 


G(s) = 


2K H cos 6 
mx x 


s 


I V I S s 2 + I V D s + 4H 2 cos 2 6 
x xx xx x 


(38) 


GH( s) = 


s 


2K K_ H cos 6 
mx fb x 


I V I S s 2 + I V D s+ 4H 2 

X XX X X X 


cos 2 6 


i 


(39) 


This system can be described, using Houpis' and D'Azzo's definitions [11], as a 
type 1 system with a step error coefficient of °°, a ramp error coefficient of 


K K 
mx fb 

2H cos 8 

x 


and a parabolic error coefficient of 0 . Next assume 


e ( s) = o . 

c 
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min 1 1 mu 



Figure 9. Uncoupled Z-axis control system. 


The forward- and open-loop transfer functions are, respectively, 


G(s) 


I 6 s + D 

XX X 


I v I S s 2 + I V D s + 4H 2 cos 2 6 

X XX XX X 


(40) 


GH( s) 


2H K K„ cos 6 
x mx fb 


I V I g s 2 + I V D s + 4H 2 cos 2 6 

X XX XX X 


(41) 


This system is described as being type 1, with a step error coefficient 00 , 

K m K fb 

a ramp error coefficient — — ■ , and a parabolic error coefficient of 0. 

r 2H cos 6 

x 
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The dual input closed-loop transfer function is 


9 (s) 
v 


2H K cos 6 9 (s) + (I S s + D ]T . (s) 

x mx c y xx x J ext v 

I V I S s 3 + I V D s 2 + 4H 2 cos 2 6s + 2H K_ K cos 6 
x xx xx x x fb mx 


•( 42 ) 


The vehicle's steady state attitude is found by applying the final value 
theorem, 


9 (t) 
v 


ss 


lim 
s— 0 
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to equation (42) . 

For a given attitude, 0 O , and a constant disturbance torque, T 0 , the 
vehicle's steady state attitude is 


9 (t) 
v 


ss 



D T 0 

+ 

2K K H cos 6 
fb mx x 


(43) 


DERIVATION OF CONTROL EQUATIONS FOR A VEHICLE 
WITH COMPENSATED PLATFORMS 


The equations are derived for a three-axis attitude control system com- 
posed of three twin compensated CMG platforms (Fig. 10) . The compensated 
platforms are obtained by attaching the platform to the vehicle's structure by a 
mechanical mechanism composed of a paralleled spring and viscous damper. 
This passive compensation mechanism is oriented so that its effects are imposed 
only on the principal control torque; the second-order cross coupling torques 
are applied directly to the vehicle. The assumptions previously used also apply 
to the compensated platforms, with the exception that the platform's coordinate 
system does not correspond to the vehicle's coordinate system. 
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A vector schematic of the X-axis platform is shown in Figure 11. Since 
the principal control torque for this platform acts about the vehicle's Z-axis, the 
compensation network will lie in the XY plane of the vehicle, thus permitting a 
relative angular displacement /3 between the platform coordinate system and the 
vehicle coordinate system. The transformation matrix relating the coordinate 
systems is 


A 

i 

P 


cos (3 

sin p 

0 


V 

A 

j P 

= 

-sin /3 

cos P 

0 


A 

3 V 

A 

k 

P 


0 

0 

1 


£ 

V 


(44) 


The subscripts p and v on the unit vectors f, 'j, and 1c indicate platform and 
vehicle coordinates, respectively. 

The angular momentum vectors H and H are described in platform 

x i x 2 

coordinates as 



A 

A 


H = 

= H cos 6 j + H 

sin 6 k 

(45) 

X 1 

x i P x i 

P 

^ 

A 

A 


H = 

-H cos 6 j +H 

sin 6 k 

(46) 

x 2 

x 2 P x 2 

P 


where <5 is the gimbal angle. 

The angular velocity of the X platform with respect to an inertial 
reference is now given as 

w * = r * j + p * j + q* k 
x X J p *x J p p 

where r *, p *, and a * are the platform components. 

X X X 


(47) 
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( 49 ) 


T = H (o * sin 6 + q * cos 6Yi - H ( r * sin 6-6 sin 
x 2 x 2 \ x tc / p x 2 \ x / 


+ H 


I (r * cos 6-6 cos 6^ 

x 2 \ x ) 


A 

k 


Summing the inverse of the torque components of equations (48) and (49) 
as was done previously and remembering that 


H 


H 

_ 

H 

Xl 


X 2 


X 


results in 


T P = T P - T P = 2H q * cos 6 
XX XjX x 2 x x x 


(50) 


T P = T P + T p = 2H r * sin 6 
xy xjy x& x x 


(51) 


T P = T P + T P = - 2H 6 cos 6 
XZ X]Z x 2 z x 


(52) 


The equation of motion for the gimbal shaft of the X-axis platform, 
obtained in the same manner as previously given, is 


T - D 6 + 2H q * cos 6 = I S 6 

mx x x be xx 


(53) 


Since the reaction torques of the torque motor and eddy current damper 
act on the platform in the 1 direction, the control torques are 


T p = D 6 - T 
xx x mx 


(54) 


T P = 2H r * sin 6 
xy xx 


(55)' 
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XZ X 

The superscript p denotes the torque components in the platform coordinate 
system. 

It is now necessary to transform these control torques into vehicle 
coordinates. This is accomplished by using the inverse of the matrix given in 
equation (44) . 


T V 


cos /3 

-sin /3 

0 


D6-T 

XX 




x mx 

T V 

_ 

sin /3 

cos (3 

0 


2H r * sin 6 

xy 




X X 

T V 


0 

0 

1 


-2H 6 cos 6 

XZ 






X 


(57) 


The control torques described by equation (57) may be simplified by 
substituting an expression that relates the vehicle rates (rpq) and the relative 
platform rate /3 to the platform rates (r *, p *, and q *) • This relationship 

X X X 

can be described as 


W x ClJ veh/ ref + 00 plat/ veh 

or in matrix form 


’r *' 


cos (3 

sin (3 

0 

X 



’x* 

= 

-sin /3 

cos /3 

0 

>* 


0 

0 

1 

m 
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Therefore, substituting these relationships into equation (57) and expanding 
yields 


T V = (d 6 
xx y x 


T V = (d 6 
xy \ x 


T \ cos £ - 2H r cos (3 sin /3 sin 6 - 2H p sin 2 /3 sin 6 
mx / x x 

(59) 


T ) sin /3 + 2H p cos j3 sin /3 sin 6 + 2H r cos 2 /3 sin 6 
mx/ x x 

(60) 


T V = -2H 6 cos 6 (61) 

xz x 


These torques are shown superimposed on the X-axis platform of Figure 12. 

Vector diagrams of the Y- and Z-axis platforms are shown in Figures 
13 and 14. The derivation of their control torques follows in a similar manner 
but with different platform to vehicle transformation matrices. The angular 
displacements of the Y and Z platforms with respect to the vehicle are ip and 
A, respectively. 

The control torques for each of these platforms in vehicle coordinates 
are as follows: 


T V = -2H /i cos M (62) 

yx y 


T V = -2H p sin fx sin a cos ol - 2H q sin 2 a sin ju + (D u - T cos a 

yy y y V y my/ 


T V = 2Hpsinusina cos 2 ol + 2H q sin a cos a sin \x 

yz r y 

+ (D ju - T 1 sin a 

\ y myy 

where ix is the gimbal angle; 


(63) 


(64) 
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Figure 12 . Schematic of the control torques generated by each compensated 

gyro platform. 


Z-axis platform 

T V = (Dy-T )sin X + 2H q sin y cos 2 A + 2H r sin y sin X cos X 
zx V z mz ) z z 


(65) 


T = - 2H y cos y 

rr 1 1 


( 66 ) 





The equations of motion for the vehicle are derived by treating each 
vehicle axis as having two degrees of freedom. 


For the vehicle's Z-axis, the equations are obtained by considering the 
X platform as one rigid body, since it is isolated from the vehicle about its 


Z-axis by the paralleled spring 


. ■ 


and viscous damper, C 


/ N-m \ 
yrad/ sec / ’ 


yrad/ sec 

and the vehicle plus the Y and Z platforms as the second rigid body. Summing 
torques about the Z-axis of the X platform and equating them to Euler's equa- 
tion of motion yields 


2H 6 cos 6 - K 

X X 


+ r * p * 


[/V - - C x[ q x* - q ] = V 

fl P - I P 1 

[_xy xx J 


(70) 


From the matrix equation (58) it can be seen that the relative platform rate is 
given as 


4 = q^ - q 

and from which the relative platform position is 

13 = f q x* ~ f q 

Substituting these expressions into equation (70) yields 

-2H 6 cos6-K/3-C/3 = I P q*+r*p*[l P -I P ~| 
x x x xz tc x x L xy xxj 


(71) 


The second equation is obtained by summing torques about the vehicle's 
axis, including the torques from the X and Z platforms as shown in 

Figure 12 and equating them to Euler's equation of motion, 
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+ K0 + C/8 + T , -2Hq sm y sin A cos A - 2H r sur A sin y 
x x ext z z 

z 

+ 2H p sin u cos 2 a + 2H q sin a cos a sin p + (D /i - T 1 sin a 

y y \ y my/ 

+ /D M - T \cos a + T' =I V p+rq|I V -I V ] 

\ y m yj ext y H H \x z ) 


( 72 ) 


The two equations for the vehicle's Y-axis are obtained by considering 
the Z platform as one rigid body and the vehicle plus X and Y platforms as 
the second. Summation of torques yields 


-2H y cos y-KA-C A = I P p* + r*q*(l P -I P ) (73) 

y r ' z z zy z z z \ zx zz / 


and 


+ K A + CA-2H p sin p sin a cos a - 2H q sin p sim a 
z z y* y 

+ 2H r sin 6 cos 2 j3 + 2H p sin 5 sin fi cos + (D 5 - T ] sin /S 
x x \ x mx j 


(v - 


T ^cos ck + T A = I ' p + rq 

y 


my/' ext_ y 




(74) 


The two equations for the vehicle's X-axis are obtained by considering 
the Y platform as one rigid body and the vehicle plus the X and Z platforms 
as the second. Summation of torques yields 


-2H /icosju-K a - C a=I P r * + p * 

y y y yx y 


a*(l P -I pN j 

y y \ y z yy / 


(75) 


and 
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+ Ka + Ca-2H r cos /3 sin /3 sin 6 - 2H p sin 2 j3 sin 6 
y y X X 


+ 2H q sin y cos A + 2H r sin y sin A cos A + (D y - T ) sin A 
z z \ z mz 1 


+ (D 6 - T Icos /3 + T =I V r+pq(I V -I V 
1 x mx/ ext x V z y 

X ' 


( 76 ) 


A block diagram of the uncoupled three-axis attitude control system is 
shown in Figure 15. This system control law, like the rigid system, employs 
only vehicle attitude feedback. 

This compensated control system can also be simplified to yield three 
sets of uncoupled equations of motion for the vehicle. Again, assume that the 
vehicle is initially at rest and that an attitude maneuver is commanded about one 
of the vehicle's principal axes. The developed reaction torques from the torque 
motor and eddy current damper associated with the controlling platform are so 
small that the vehicle's velocity about that axis essentially remains zero. Also, 
the angular displacement between the vehicle and platform is very small. Hence 
equations (70) through (76) are simplified to the following: 


X-axis 


- 2H /i cos n - K a - C o = I P r * 

y y y yx y 


+ Ka + Ca+T . =I V r 
y y ext x 

y 


Y-axis 


-2H y cos y - K A-C A = I P p * 
z z z zy z 


+ KA + CA + T ,=I V p 
z z ext y 

t r * 


(77) 

(78) 


(79) 

(80) 
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PLATFORM EQUATIONS OF MOTION y 
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VEHICLES EQUATIONS OF MOTIONS 


VEHICLE RATE 


P' T x"* K y ai C y a-2H x r sin 8 sin ft cos /3 - 2H X psin 8sin 2 /3+2H 2 qsiny cos 2 X+2H z r sinysin Xcos X+(D Z y-T m )sin X+(D x 8-T m )cos/3+T #xt -pqQj-lX )] X 

2/rJ-+K z X + C z X-2H y psin^. sin a cos a - 2H y q sin/i sin 2 a+2H x r sin 8 cos 2 /?+2H x psin 8 sin /3co$ /3+(D x S-T m )sin/3+(D y /i-T m )cos a +T ext -rq(I x -I* ) ! J_ _p 
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Figure 15. Three-axis compensated platform control system. 












Z-axis 


2H 6 cos 6 - K j 8 - C j3 = I P q* 

xz 


X 


X 


+ K/3 + C/3+T = I V q 

x x ext z 

z 


( 81 ) 

(82) 


It is now possible to discuss and analyze control for each axis independ- 
ently. Figure 16 is a block diagram of the uncoupled Z-axis control system. 
This, like the rigid system, is also a dual input control system. Assuming that 

T = 0, the forward- and open-loop transfer functions are, respectively, 
ext 


G(s) = 


2K H cos 6|C s + K 1 
mx x Lx xj 



2K H ( 


mx x 

S 1 

I V I S I P s' 

l 

X XX xz 

I v 

I P D 1 s 3 + 

X 

xz xj 


fl g C (l r +l e ) 

XX X\ X xz / 


I S K +C D U 4H * I v cos 2 5 1 s 


XX X XX 




2 t V 


X X 


K D (l v + I pN \ 
X x\ X xz J 


I V + I P 1+4H 2 C cos 2 6 Is + 4H 2 K cos 2 6 


x xz I xx 


5 £ 


X X 


GH( s) = 


2K K H cosfifc S+K ) 
mx fb x \ x x/ 


s-|I V I S I 

X XX XZ 


p g C (i + I P Vl V I P d 1 a 

:z I xx x\x xzy x xz xj 

f(l V + I pN \ (l S K +C D \+ 4H. 2 lV cos 2 6l s 2 
[Vx xz M xx x xxi xx I 


K D (i + I P 1 + 4H 2 C cos 2 6 s + 4H 2 K cos 2 6 

xxyxxz^ XX XX 


4 


(83) 


(84) 


This system is described as a type 1 with a step error coefficient of °° , 

K mx K fb 

a ramp error coefficient of —rr , and a parabolic error coefficient of 0 . 

2H cos 6 
x 
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Assuming that 6 (s) = 0, the forward- and open-loop transfer functions 
are, respectively, 


St P o3,/t S 


I 6 I F s 3 +(l g C +1 P D V + ( ] 

XX XZ \ XX X XZ X/ \ 


I & K + D C 

XX X XX 


G(s) = 


+ 4H 3 cos 2 S )s + K D 

X / X > 


s ;i v i s i E v + 

X XX xz 


fl S C (l V + I P V I v I P D 1 s 

I xx x\ x xz / x xz xj 
\(l'+l*Vl*K +c D \+ 4H 2 I v cos 2 $ls 2 

|_\ X XZ/V XX X x x/ XX J 


[ K x 


I P I S C s 4 + 

XZ XX X 


I V + I P ^+4H 2 C cos 2 6 s + 4H 2 K cos 2 6 f (85) 
xxz/ xx J xx ‘ 


/l P D 

y xz 

$js 2 + ^ 


C +I P I S k'\s 3 + 

XX XZ XX X 


/l P D 

y xz 


K 
x x 


H cos 6 ) 
: x ) 


+ 4H 2 C cos 2 6 Is 2 + (4H 2 K cos 2 $ + 2K_K C H cos 0 I s 
xx r \ x x fbmxx 


+ 2Kj_ K H K cos 6 
GH( s) = — . » x x 

s I V I S I P s‘ + ri g C fl V + I P V I v I P D 

X XX XZ I XX x\ X XZ/ X XZ X 


‘v4i g c fi v + i p > 

:z xx x\ x xz / 

[(i v + I P Vi S K + C D \+ 4H 2 I v cos 2 6 
l\x xz/yxx x x x/ xx 

[k D (l V + I P \ + 4H 2 C cos 2 6 s + 4H 2 K cos 2 6 
x xyx xz/ xx J xx 


(86) 


This system is described as a type 1, with a step error coefficient of 


K K 
mx fb 


°°, a ramp error coefficient , and a parabolic error coefficient of 0. 

’ 2H cos 5 
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The dual input closed -loop transfer function is 


2H K cos6(C s +K )0 (s) +|I g I P s 3 +/l g C 
x mx y x x J c I xx xz l xx x 

+ I P D )s 2 +(l S K +D C + 4H 2 cos 2 6 )s +K D T (s) 
q / x_ xz x/ \xx x _x x x / x xj ext 

vS = I v I p I*s s + l7l v + I P '\l g C + I v I p D 1 s 4 
x xz xx y x xz I xx x x xz xl 


‘J 


+ 

Yi v + 1 pN ) 

(l g K + C D N 

1+ 4H 2 I V cos 2 8 

s 3 + 


\ X XZ J 

y XX X X x / 

1 XX 



X 


+ 4H 2 C cos 2 6 s 2 + 2H K_ K C cos 8 + 4H 2 K cos 2 $ I s 
xx I I x fb mx x xx 1 


+ 2H K K _ K cos 6 
x mx fb x 


(87) 


The vehicle’s steady state attitude for a given attitude command, 0 O , 
and a constant disturbance torque, T 0 , is 


D T n 


6 (t) = 

v ss K 


fb 


2K K_ H cos 6 
mx fb x 


(88) 


SINGLE-AXIS ANALYSIS 


The advantages and disadvantages of the compensated platform system 
with respect to the rigid platform system are determined by analytical methods. 
This analysis was performed principally by use of an analog computer; however, 
a root locus analysis was also performed on the linearized systems. 

Since the analytical results obtained here are to be verified by experi- 
mental simulation of the control system, the "physical" parameters used in 
the analysis are those of the simulator: 
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I V = 271 kg-m 2 
I P = 0.271 kg-m 2 
I S = 0.001355 kg-m 2 
H = 0.676 N-m-s 

The other "performance” parameters (K, C, D, K , and K_ ) were determined 

m id 

so that the system’s transient response for a one -degree step command cor- 
responds to that of a critically damped system. This was accomplished by 
programing the single-axis compensated platform control system, shown in 
Figure 16, on an Applied Dynamics 8800 analog computer. 

This computer is equipped with a "rep-op" mode of operation. This mode 
time scales the problem so that its transient response occurs in a very short 
time. The computer is then adjusted to switch automatically from reset to oper- 
ate at a frequency that corresponds to the settling time of the time-scaled prob- 
lem. When the system's transient response is viewed on an oscilloscope, it 
appears as a continuous trace. Thus, watching the trace and varying the poten- 
tiometer, which corresponds to the respective "performance" parameters, gives 
the desired transient response. 

In addition to satisfying the imposed system performance requirement, 
several other conditions must also be fulfilled: 

1. The gains (K^, K^) must be such that the gimbal angles do not 

exceed or become saturated at either ±60 degrees, during the one-degree 
maneuver. 

2. The magnitude of the spring rate (K) and viscous damping coef- 
ficients (C,D) must be realistic to enable implementation on the simulator. 

3. To minimize the attitude error induced by external disturbance 

D 

torques, the term — — — — — — should be kept as small as possible. 

2H K K_ cos 6 
x my fb 

This can be seen from equation (88) . 
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With these constraints in mind, the following "performance" parameters 
were obtained: 

K = 0.271 N-m/rad 
m 

K_ =0.8 rad/rad 
fb 

K = 135. 5 N-m/rad 

„„ _ N-m 

C — 67.7 , / 

rad/ sec 


D 


0.014 


N-m 
rad/ sec 


The results of the analog studies for the compensated platform control 
system are shown in Figure 17. The system's performance curves confirming 
the specified transient response are seen in Figure 17 ( Sheet 1) . The time 
constant and settling time were measured from the vehicle position trace to be 
7.00 seconds and 25.80 seconds, respectively. The peak torque developed during 
this maneuver was 0. 189 N-m. This system's response to a step command of 
0. 04 radian is seen in Figure 17 (Sheet 2) . This is the maximum attitude maneu- 
ver this system can perform without exceeding or saturating the ±60-degree 
gimbal angle constraint. 

The electrical power required to drive the gimbal shaft torque motor was 
calculated using the following equation: 


Power = 


K_ 6 
fb v 


(89) 


The constant G is determined from the electrical characteristic of the torque 

(volts/rad) 2 


motor. It has the units 


ohm 


The power curve, shown in Figure 18, 


was derived for the one-degree maneuver by using a value of G= 454 and 
measuring 0 from the vehicle trace in Figure 17 (Sheet 1) . This particular 

value of G was used because it corresponds to the torquer incorporated in the 
experimental control moment gyro. Thus the calculated power curve can be 
compared with the recorded experimental power curve. 
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The energy consumed by the torquer was calculated using the following 
equation . 


N ' / 
Energy = £ P ave At. = J 

i= 1 i 0 


Power dt 


(90) 


The energy curve shown in Figure 19 was derived by measuring the area under 
its power curve, Figure 18, with a K&E planiometer. 

A root locus analysis was also performed on this system. But, because 
of the variable term cos 6, which appears in the open-loop transfer function, 
the analysis is valid only for the particular gimbal angle for which the equation 
is solved. Choosing the gimbal angle as zero and substituting the values of the 
parameters found from the analog study into equation ( 84) yields 

K ( s + 2) 

GH(S) " s(s 4 +’261s 5 + 8075s r +~6425s + 2500) (91 


where K_ is defined as the static loop sensitivity. The roots of the denominator 
were found by using an IBM 94 digital computer. They appear in the open-loop 
transfer function as 


K_( s + 2) 

GH ^ S) = s (s + 35To8)~ (s+ 224.72) ( s + 0. 402£T+~To. 39347 Ts + 0.4029 - j 0.3933) 

(92) 

The effects on the predominant open-loop poles and zero caused by various 
spring rates (K) and viscous damping coefficients (C) are shown in Figure 20 
for a gimbal angle of zero degree . 

The root locus of equation (92) is shown in Figure 21. The closed-loop 
poles correspond to a static loop sensitivity of 200, which was determined by the 
analog study. Also seen in Figure 21 is the effect of the cos <5 term on the predom- 
inant (complexed) closed -loop poles, ijhe damping ratio and natural frequency 
of these poles vary from 0.275 < £ < 0.063 and 0.2 < w^(rad/sec) < 0.488, 

respectively, for gimbal angle excursions of 0 to 60 degrees. 
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CO 


Figure 20. Open-loop poles and zero values for various spring and damping coefficients. 





The closed-loop transfer function of this dual input system, for a gimbal 
angle of zero degrees, appears in factored form as 


e (s) 

V 


25O(s+2)0 (s)Mf.005(s + 0.65) (s+ 35.29) (s+ 224.4) (s) 

( s + o. 25) ( s + 35.' 08) ( s + 224. 72) ( s 2 + 0. 555s + 0. 2012) ' 

(93) 


Bode plots of 


Mjw) 

0 (jw) 

c 


and 


Mjw) 


are shown in Figures 22 and 23, respec- 


tively. The frequency bandwidth of this system is 0.25 rad/ sec. 


The equation of motion for the vehicle in the time domain, for a step 
command A 0 , assuming the external disturbance torque is zero, is 


0 y (t)= A 0 [l. 25 - 1. 79e °’ 25t _ 29.4 x 10 6 e 35 - 08t 

+ 1.08 x 10~ 6 e ~ 224 * 72t _ 1.0l8e“° -276t sin(0.3535t - 31.8°)] 


(94) 


The transient response curves of the rigid platform control system are 
used to form the basis upon which the effects of the compensation network can be 
determined. For this comparison to be valid, the control system shown in 
Figure 9 was also programed on the Applied Dynamics 8800 analog computer and 
the same system parameters and attitude commands were applied to this system 
as to the compensated systems. 


This system's transient response curves for a step command of one 
degree are shown by the solid lines in Figure 24. The measured time constant 
and settling time for this system are 7. 15 and 19.40 seconds, respectively. The 
peak developed control torque was 0. 298 N-m. The dotted lines superimposed in 
this figure are the response curves of the compensated platform system for a 
step command of one degree. Thus, the difference between the performance 
characteristics of the two platform configurations is easily distinguished. These 
curves were traced from the original strip chart recordings. 


The power required to drive the gimbal shaft torque motor during this 
one-degree maneuver was also calculated using equation ( 85) . The value of G 
was 454 and 6 was again measured from the vehicle’s attitude trace in Figure 

24. This system's power curve is also shown in Figure 18. The corresponding 
energy curve was developed by measuring the area under its power curve, which 
is shown in Figure 19. 
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A root locus analysis was also performed on this system. Like the 
compensated system, this system also has the variable term cos <5 in its open- 
loop transfer function. Choosing the gimbal angle to be zero and substituting 
the values of the parameters into equation (36) yields the open-loop transfer 
function 


K 

GH( S) ~s‘( s + 0. 5074) ( s + 9. 852) 


(95) 


The root locus of this system is shown in Figure 25. The closed-loop 
poles correspond to a static loop sensitivity of 0.8, which was defined by the 
analog study. The effects of the cos 8 term on the complex poles are also shown 
in Figure 25. These poles have a damping ratio and natural frequency variance 
of 0.3 < £ < 0.878 and 0. 197 < ^(rad/sec) < 0.285, respectively, for gimbal 

angle positions of 0 to 60 degrees. 

The closed-loop transfer function in factored form for this dual input 
system for which the gimbal angle is 0 degrees is 


9 (t) 
v 


6 ( s) + 0. 005 ( s + 10. 36) T ( s) 

C GXt 

(s + 9. 86) (s 2 + 0. 498 s+ 0. 0809) 


(96) 


Assuming that the external disturbance torque is zero, the equation of 
motion of the vehicle in the time domain for a step command A 0 is 


Mt) = A 0 £l. 25 - 1. 098 x 10" 3 e~ 9 ° 86t 

sin( 0. 1375t + 29. 4°)J 


- 2. 59e 


-0. 245t 


(97) 


Bode plots of 


0 v (j") 

l V5“> 


and 


® v (3^) 




are shown in Figures 26 and 27, 


respectively. This system has a frequency bandwidth of 0.245 rad/ sec. 
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Figure 25. Root locus of rigid platform system. 


Figure 26. Normalized Bode plot 
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Figure 27. Normalized Bode plot for the rigid platform system. 



NONLINEAR CONTROL LAW 


During the analog study previously performed, it was apparent that the 
gimbal angle constraint of ±60 degrees without saturation seriously restricted the 
magnitude of the attitude maneuver. With this fact in mind, a nonlinear control 
law was developed that enables the same control systems to perform attitude 
maneuvers of up to 360 degrees at rates exceeding those developed by the linear 
systems previously analyzed. The control system's performance will be demon- 
strated by analog and experimental simulation. All of the previous assumptions 
and constraints are still valid with the exception of the gimbal angle saturation 
constraint. 

This control law produces performance characteristics similar to those 
of a CMG-reaction control jet system; that is, for attitude errors exceeding the 
chosen limits, ±A0 O , the generated control torque approximates an impulse. As 
the error is decreased to within the ±A0 O limits, the control torque becomes 
continuous and proportional, thereby enabling the system to acquire and main- 
tain any given attitude command. 

A block diagram of the nonlinear control law and the gyro-recaging 
(momentum dump) logic is shown in Figure 28. An explanation and analog 
simulation of the recaging logic will be presented later. An explanation of the 
nonlinear control law follows. 

Assume that the vehicle's initial conditions and external disturbance 
torques are zero and that a positive step command, A^, has been applied to the 
system. The developed attitude error, -A0(t) , is applied to a nonlinear func- 
tion, which is described as having the following piecewise characteristics: 


v(t) = A0( t) 

-A0 O < A0(t) A0 O 

v( t) = V 

o 

A0(t) > AO 0 

v(t) =-V 

o 

A6( t) < - A0 O 


Assuming that the attitude error is greater than -A0 O , the output of the nonlin 
function is -V volts. This voltage is in turn applied to a relay that is driven 

° I. 

by the output of logic block III. This relay is polarized so that it is closed for a 


ear 
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Figure 28. Nonlinear control law. 










logic level input of i and open for 0. At this instant the attitude error is negative 

and the gimbal angle is less than -60 degrees. Thus the outputs of logic blocks 

I and II are -1 and 0, respectively. Their sum (-1) is applied to logic block III, 

thus producing an output of 1. Therefore the relay is closed and the voltage 

(-V ) is applied to the torque motor. Since this voltage is greater than that 
o 

which would be applied during linear control, the related motor torque, gimbal 
rate, control torque, and vehicle rate are also greater. When the gimbal angle 
reaches -60 degrees, the brake logic, also shown in Figure 28, induces an 
electromagnetic brake located on one of the gimbal shafts to be energized, there- 
by holding the gimbal shaft at -60 degrees. At the same time the output of logic 
block II changes from 0 to -2. The sum, now -3, produces an output of 0 from 
logic block III. This opens the relay and removes the voltage (-V^) from the 

torque motor. The vehicle continues to rotate at its developed angular velocity 
until the vehicle's attitude exceeds the commanded attitude. At this time the 
brake is de-energized and the output of logic block I changes from -1 to +1. The 
sum (-1) produces an output of i from logic block III. Thus the relay is again 
closed and the error signal, now within the ±A 0 limits, is applied to the torque 
motor. The control torque is now continuous and proportional, thus enabling the 
system to acquire and maintain the commanded A 0 attitude. 

This nonlinear control law was applied to both the rigid and the compen- 
sated platform systems. The performance of these systems was simulated using 
the Applied Dynamics 8800 analog computer. The piecewise characteristics of 
the nonlinear function were arbitrarily chosen to be 

v(t) = A0(t) - 0. 0175 rad ^ A0(t) ^ 0. 0175 rad 

v(t) = 0.2 volts A0(t) < 0.0175 rad 

v(t) =-0. 2 volts A0( t) <-0. 0175 rad 

while the system's parameters correspond to those used in the analog studies. 

The performance curves for the systems employing the nonlinear control 
law are identical to their corresponding system configurations studied for 
attitude errors equal to or less than 0.0175 rad. 

The transient response of the compensated platform system to step 
commands of 0. 175 rad is seen in Figure 29. The peak control torque developed 
during this maneuver was 2. 57 N-m. 
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The electrical power applied to the gimbal shaft torque motor was cal- 
culated by using the following equations during nonlinear time of operation 

Power = V /R 2 (98) 

o 

where R is the motor's armature resistance and by using equation (89) during 
linear operation. A value of G = 454 and R = 12 ohms was used to derive the 
power curve in Figure 30, since the same control moment gyro will be used to 
perform all the various system and control law experimental maneuvers. The 
attitude of the vehicle was read from the vehicle's attitude trace in Figure 29. 

The energy curve shown in Figure 31 was derived by using equation (90) 
and its power curve in Figure 30. The transient response curves for a step 
command of 0. 175 rad applied to the rigid platform control system, which 
employs the same nonlinear control law, are shown by the solid lines in Figure 
32. The dotted lines indicate the response of the compensated platform, non- 
linear control law system shown in Figure 29. The difference between the per- 
formance characteristics for each system is clearly visible. These curves were 
also traced from the original strip chart recordings. 

The torque motor's power and energy curves were derived by using the 
same equations and constants that were used for the compensated platform 
system. The vehicle's attitude was measured in Figure 32. These curves are 
also shown in Figures 30 and 3 1 . 

A rather unique feature of this nonlinear control law, as seen from the 
analog simulation of both systems, is that for attitude maneuvers of greater than 
one degree, the power and energy requirements are a constant; that is, they are 
independent of the magnitude of the maneuver. 

To explain the function of the recage logic shown in Figure 28, it is 
sometimes easier to think of the control moment gyro as a momentum manage- 
ment device; that is, it can absorb either positive or negative angular momentum 
from the vehicle. The amount of manageable angular momentum is a function of 
the magnitude of the gyro's angular momentum and gimbal angle constraint. For 
the systems under consideration, it is found to be 

H ,, = 2Hsin 6 = (2) (0.676) sin 60° = 1. 17 N-m-s 

manageable max 
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The presence of an external disturbance torque on the vehicle induces 
angular momentum into the system. This additional angular momentum must be 
absorbed (cancelled) if the angular momentum of the system is to remain con- 
stant and thus maintain the vehicle in a nonrotating state. This is accomplished 
by precession of the gyro's angular momentum vectors in the direction that will 
nullify the added angular momentum. This process appears in equation form as 

T t + 2Hsin 5(t) = 0 . (99) 

GXt 


The angular velocity of the vehicle is assumed to be zero. 

When the gimbal angle reaches its designated constraint, the gyros are 
no longer capable of absorbing any additional angular momentum. The recage 
logic is, therefore, devised to remove this absorbed angular momentum so that 
the gyros can again perform their momentum management tasks. This is nor- 
mally accomplished by firing a set of reaction control jets. 


The recaging logic is as follows: When the gimbal angle becomes satu- 
rated at either ±60 degrees and the rate of change of the absolute value of the 

attitude error ^ is positive, the firing command is given to the appropriate 

ut 

pair of reaction control jets. 

The length of time the thrusters are fired is determined by the following 
equation 


2H sin 6 

max 


2T t 


( 100 ) 


where the term T is the torque developed by a control jet and 6 is the 
t max 

gimbal angle constraint. 

An analog simulation of the recaging process for the rigid and compensated 
platform systems is shown in Figures 33 and 34, respectively. The disturbance 
torque for both cases was 6. 78 x 10 N-m. The 1. 17 N-m-s of stored angular 
momentum was removed by firing 0.469 N-m thrusters for 2.5 seconds. 
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Figure 33. Rigid platform — = 0.00678 N-m. 


EXPERIMENTAL SINGLE-AXIS SIMULATION 


A single-axis experimental simulator was built to verify the analog data 
previously presented. This simulator (Fig. 35) was designed with the capabilities 
of performing maximum attitude maneuvers of ±10 degrees, with an accuracy 
of 0. 1 degree. The simulation hardware will be discussed with respect to three 
areas: the twin CMG, the single-axis vehicle motion simulator, and the control 
and recording electronics . 


61 


I 






Figure 35 . Single-axis two-degree-of-freedom motion simulator. 

Twin Control Moment Gyro 

The design of the twin CMG was influenced primarily by two factors: the 
first is the physical parameter; that is, the gimbal shaft’s mass moments of 
inertia, the gyro platform’s mass inertia, the gimbal shaft damping coefficient, 
and the angular momentum of the gyros; and the second is the desire to maintain 
the lowest possible friction level about the gimbal shafts. This desire for 
minimal gimbal shaft friction is based on the knowledge [5] that pointing accu- 
racy (attitude error) is inversely proportional to the gimbal shaft friction. Fig- 
ure 36 is a photograph of the experimental twin CMG. 
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The gyro rotors and housing were obtained from two Minneapolis Honey- 
well JG044D-4 cageable vertical gyros. The gyro motors were rated as having 
0.630 N-m-s angular momentum. The spin motors are two phase 400-Hz, 
115-volt synchronous motors, which rotate at 21 000 rpm. The power to drive 
these spin motors was obtained from two number 42 wires that were looped down 
from the wire support bracket and passed through the hollowed gimbal shafts. 

A 0 . 33 /iF capacitor was cemented on the housing of each gyro to develop the 
two-phase voltage required by the spin motors. The required voltage was there- 
by obtained from two wires instead of the normal three; thus the imposed wire 
spring rate acting on the gimbal shafts was reduced by 33 percent. 

The gimbal shaft torque motor was a brushless dc torquer, Model 
#TQ18W-23, manufactured by the Aero Flex Corporation. This motor features 
no commutator, brushes, or contact friction. It has infinite resolution; that is, 
a linear torque output versus current input over a range of ±60 degrees. This 
motor has a permanent magnet rotor and will develop a peak torque of 0. 845 N-m. 
Its electrical characteristics are a motor sensitivity of 0. 705 N-m/amp, an 
electrical time constant (L/R) of 8 x 10 4 sec , and a developed back emf of 
0 . 07 volts/ rad/ sec . 

An eddy current damper was incorporated into the CMG to develop the 
needed gimbal rate damping necessary for stability. This damper was designed 
and built by the Astrionics Laboratory of Marshall Space Flight Center. The 
damping coefficient was developed by rotating a 0. 208 cm copper disk between 
eight electromagnetic poles, which were equally spaced around an 8.14-cm- 
diameter circle. After fabrication of the unit, calibration tests dictated that 

0. 86 amps were required to develop the 0.014 — — damping coefficient. 

3? SCI/ S6C 

Its corresponding power consumption was 9.26 watts. 

The gimbal shaft potentiometer was manufactured by the Bendix Cor- 
poration. This device operates on the reluctance principle and therefore has no 
windings on the rotor. Thus the necessity for slip rings or brushes was elimi- 3 
nated and thereby reduced the unit's friction to the breakaway level (4.53 x 10 
N-cm) of its bearings. This unit requires 10-volt, 400-Hz power and yields a 
linear output over a range of ±60 degrees with an accuracy of less than 2 degrees 
error. However, an accuracy of 20 arc min is obtainable over the range of 
±45 degrees. 
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The gimbal shaft brake is a Model BA- 100 made by Dial Products 
Company. This is an electromagnetic brake designed for on-off service. The 
brake is composed of two separate pieces, the electromagnet and the pole shoe 
assembly. The pole shoe assembly consists of a piece of magnetic iron attached 
to the gimbal shaft hub by a spring disk of heat-treated beryllium copper. This 
unit features zero backlash, zero residual drag, up to 3 degrees allowable 
angular misalignment, and response time of approximately 10 ms. The brake 
is energized by 28 mA and 100 volts dc. 

The two aluminum spur gears used to connect the two gimbal shafts were 
machined by the Pic Design Corporation. These are precision machined gears 
with 288 teeth, a pressure angle of 20 degrees, and a pitch diameter of 10. 16cm. 
During assembly of the twin CMG, the gears were mated so that minimum back- 
lash (approximately 30 arc sec) and contact friction were realized. 

The gimbal shaft bearings were manufactured by New Hampshire Ball 
Bearing, Inc. These are instrument-type ball bearings and were machined to 
ABEC-7 tolerances. These bearings are composed of 8-400C stainless steel 
balls held by a stainless steel crown-type retaining ring. No side shields or 
lubrication was used with these bearings. A nominal breakaway friction level for 
these bearings is 2. 96 x 10 N-cm. They were mounted in their housings so that 
one bearing on each shaft was permitted to float while the other was held fixed. 
This configuration eliminates the axial loads that would occur as a result of 
thermal expansions. 

Another feature incorporated into the gimbal shaft design was the use of 
two flex couplings for minimizing bearing side loads that would result from 
misalignment of the gimbal shafts . The machining tolerance on all dimensions 
related to alignment and bearing fits were held to ±0.005 mm. All of the above 
precautions were in keeping with the design objective for minimum gimbal shaft 
friction. 

The electrical power was supplied to all of the CMG components by a 
conical arrangement, consisting of 22 number 41 copper-varnish coated wires. 
These wires were suspended from a circular ring attached to the ceiling. This 
method of power transmission enables the simulator to perform attitude maneu- 
vers of ±10 degrees with no measurable loading effects. 
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Single-Axis Motion Simulator 

The experimental simulation of the vehicle performing single-axis 
maneuvers was accomplished primarily through the use of a single -axis, two- 
degree-of-freedom, air bearing motion simulator. The test capabilities and 
design objectives to be met by this simulator were as follows: 

1. It must have the same physical characteristics previously given; 
that is, mass moment of inertia of the vehicle and the gyro platform. 

2. It must be as free of contact friction and external disturbance torques 
as is practical. 

3. It must be capable of simulating either the rigid or the compensated 
platform configuration. 

4. It must be capable of performing momentum dumping operations. 

The first design objective was fulfilled for the gyro platform by sizing 
(weight) and locating the various CMG components on the gyro platform. In 
addition, the center of gravity of the gyro platform was made to lie on the sim- 
ulator's axis of rotation. To satisfy the mass moment of inertia requirement 
for the vehicle, a structure consisting of four cantilever arms was built, and 
to each arm was attached a 38. 5 kg mass. The center of gravity of the simulated 
vehicle was also made to lie on the axis of rotation by placing small lead weights 
on the cantilever arms. 

The second and third items were accomplished by the use of two concentric, 
combination "thrust -journal" air bearings shown in Figures 37 and 38. The 
outer- or vehicle -bearing supports the simulated vehicle inertia and also isolates 
it from the simulator stand. This bearing is composed of both a thrust bearing 
and a journal bearing. The thrust bearing has 16 nozzles equally spaced around 
a 21.05 cm-diameter circle. Each nozzle has an orifice diameter of 0.356 mm 
and is pressurized by a common manifold to 2.06 x 10 5 N/m 2 gage. The journal 
bearing was incorporated into the design to produce a high torsional stiffness in 
the plane perpendicular to the journal bearing's axis of rotation. This is a 
15.20-cm diameter bearing machined to have a 0.0508 ± 0.00508-mm diametrical 
clearance. The bearing is composed of two rows, each having 12 nozzles equally 
spaced around the circumference. These nozzles were designed specifically to 
develop a high stiffness ( 12. 6 x 10 6 N/m) . This was accomplished by using an 
orifice diameter of 0.305 m and pressurizing them to 5.05 x 10 5 N/m 2 . A set of 
three equally spaced vent holes was machined into the bearing at the intersection 




Figure 37. Two -degree -of -freedom single-axis air bearing simulator 




of the thrust and journal bearing surfaces. This was done to stabilize the bearing 
by eliminating any cross coupling of journal bearing air pressure with the thrust 
bearing air pressure and also to simplify the design analysis [ 12, 13] . 

The inner or platform bearing is used to form the compensated gyro 
platform configuration. This bearing is also a thrust-journal bearing. When 
pressurized, this bearing supports the CMG platform and isolates it from the 
vehicle. The torques generated by die twin CMG are then transmitted through the 
passive compensation network to the vehicle. To maintain a minimum contact 
friction level, a labyrinth seal was employed to supply air to the bearing. This 
dictates the use of one common source of air pressure for both the thrust and 
journal portions of this bearing. This thrust bearing was also composed of 16 
nozzles equally spaced around a 15.20-cm diameter circle. Each of these noz- 
zles has a 0. 177 -mm diameter orifice and is pressurized to 2.06 x 10 5 N/m 2 gage. 
The journal bearing has 15. 20 -cm diameter and was also machined to have a 

0. 0508 ± 0.00508-mm diametrical clearance. This bearing is also composed of 
two rows of 12 nozzles, each orifice having a diameter of 0.305 mm. These 
nozzles develop a stiffness of 2.72 N/m when pressurized to 2.06 x 10 5 N/m 2 gage. 

The air pressure for each of the bearing's supply manifolds was regulated 
by a one-stage Grove pressure regulator. The air from each of three pressure 
regulators is then filtered by a 10 -micron absolute filter before entering the 
bearings. This precaution prevents foreign particles from clogging the nozzle 
orifices or from scratching and galling the bearing surfaces. 

The passive compensation network is shown in the photograph in Figure 
36 and by a sectioned illustration in Figure 39 . It should be noted that two 
compensating systems were employed, one on either side of the gyro platform. 

This configuration tends to induce pure rotation of the platform, with respect to 
the vehicle, by cancelling unbalanced forces resulting from mechanical misalign- 
ment of the compensating devices. 

The compensation springs are helical with a linear spring rate of 

1. 07 x 10 3 N/m. The spring housing was designed to preload the springs so that 
even during their deflections they were always kept under compression. This 
method ensures a linear force versus displacement relationship. 

The viscous dampers were manufactured by Houdaille Industries, Inc. 

These are piston-type, linear motion dampers. The dampers were modified by 
removing the unit's piston shaft O-ring seals and operating them submerged in a 
bath of silicon based oil. This peculiar method of operation was necessary to 
eliminate the high coulomb friction level induced by the O-rings. Silicon oil 
was chosen as the damping fluid for its low thermal viscosity gradient. 



Figure 39. 


Passive compensation network. 



The fourth simulator design objective concerns the capability of per- 
forming momentum dumping operations. This was met through the use of an 
ac torque motor. This motor was incorporated into the design of the vehicle air 
bearing by bolting the rotor, which consists of a solid steel cup, to the bottom of 
the journal portion of the vehicle air bearing. The stator was in turn mounted 
into the base of the fixed half of the vehicle's air bearing manifold. Thus a pure 
torque was applied on the vehicle in the same manner as reaction control jets 
would be applied on an actual vehicle, and yet there are no brushes, slip rings, 
or other contact type devices to impart retarding torques that would inhibit 
normal simulation maneuvers. 

The gimbal shaft torque motor was an Inland Model X-3001. This motor 
operates on two-phase, 400-Hz power and develops a peak torque of 0.472 N/m. 
During simulation the fixed phase voltage was set for 115 volts, while the control 
phase voltage was set so that the motor would develop a torque of 2 T , which 

t 

would satisfy equation (99) . 


Electronics 

The third major area to be discussed concerns the associated electronic 
devices used to implement and record this control system. Figure 40 is a 
photograph of the equipment. 

The system gains, and K , were obtained from two Hewlett- 

Packard Model 2470A data amplifiers. These amplifiers were chosen for their 
gain stability (±0.005 percent per month) and output linearity (±0.002 percent 
of full scale) . These amplifiers have a maximum output of 100 mA ± 10 volts. 

Additional amplifiers were obtained from two analog computers; a 
Pace-TR-10, which has 20 amplifiers with a maximum output of ±10 volts, and 
a Donner-3200, which has 10 amplifiers with a maximum output of ±100 volts. 

The Pace amplifiers were wired to act as comparators. The comparators were 
set to develop a +5. 0-volt output. These outputs were used to drive digital logic 
elements. The Donner amplifiers were wired to differentiate and filter the 
vehicle, platform, and gimbal position signals to yield their respective velocities. 
These velocities were derived solely for the purpose of recording. 

The nonlinear control law's logic was developed by using Sylvania SNG 
14/F-530 NAND gates and Siemens polarized relays. The NAND gates were 
triggered by the +5 volts of the comparators, and the relays were in turn 
energized by the outputs of the NAND gates. These relays have a response time 
of 10 ms. 
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The recordings of both analog and experimental simulations were made 
on an eight -channel, Model RF-1783, Brush recorder, which has a dc linearity 
of 0.5 percent full scale, an input impedance of 1.25 megohms, and a time con- 
stant of 2 . 65 ms . 

Analog measurements of the vehicle's and gyro platform's attitudes were 
made with an electro-optical measuring system. The vehicle's attitude measure- 
ment was used as the feedback signal for the control system. Both measurements 
were also used to record their respective transient behavior during maneuvers. 

The system functioned by reflecting the light from a fixed high intensity 
light source off a mirror mounted on the side of the moving vehicle and into a 
lens. The lens converged the light to a point on the face of a photovoltaic ele- 
ment. The element in turn develops an output voltage proportional to both the 
light's intensity and position. Figure 35 shows both the vehicle's and platform's 
attitude measuring system. The light source was a Sylvania FAL High Silical 
Halogen lamp. This lamp develops 11 000 lm, which is enough to saturate the 
intensity variable of the photovoltaic element. Thus, the attitude sensor's output 
was reduced to only a function of the position of the light source's image. The 
photovoltaic device is called a Radiation Tracking Transducer (R.T.T.) and 
was manufactured by Electro-Optic Systems, Inc. This unit is rated by the 
manufacturer to have an output responsitivity of 19. 7 mV/mm, an output linearity 
of 2 percent, and a usable angular field of ±6 degrees with a standard lens. The 
Angenieux lens has a focal length of 2.54 cm and f/0.95. 

A schematic of the vehicle's attitude measuring system is shown in 
Figure 41. The following mathematical relationships were derived to establish 
the distances W and L at which the light source and R.T.T. must be placed 
from the mirror for the system to be capable of measuring vehicle angles of 
±10 degrees. 

The zero position of the measuring system is adjusted so that the R.T.T. 
and the light source lie along a line of centers which runs through the axis of 
rotation of the vehicle and the longitudinal center of the mirror. The reflecting 
surface of the mirror is adjusted to be perpendicular to this imaginary line. The 
dimension AR, shown in Figure 41, defines the distance between the mirror's 
surface and the edge of the vehicle simulator, measured along the line of centers. 
This dimension is given by the geometrical relationship, 


AR = 


R( 1 - cos 6) 
cos 9 


( 101 ) 
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Figure 41. Schematic of the attitude tracking system. 

where R is the radius of the vehicle disk. The distance x> shown in Figures 
41 and 42, is the lateral displacement of the light-source image when measured 
from the line of centers. By the use of Figure 42, two geometric relationships 
for x were obtained, 
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X = t(W - AR) cos 20] tan 2 0 


( 102 ) 


X = [(L - AR + (W - AR) cos 20] tan a . (103) 

Substituting equation ( 10 1) into equations ( 102) and ( 103) and then equating 
equation ( 102) to equation ( 103) yields 


L = 


W - R 


( 1-cos 0) 
cos 0 


tan 2 0 


tan a. 


cos 20 + R 


( 1-cos 0) 
cos 0 


(104) 


The equation relating W and L for the vehicle's attitude measuring 
system was obtained by substituting the following parameters into equation ( 104) : 

R = 38. 1 cm 

0 =10.5 deg 

a = 6.0 deg , 

therefore 

L = 2.476 W + 0. 1952 . (105) 

Equation ( 105) was satisfied by letting W equal 25.4 cm and L equal 63.4 cm. 

The dimensions W and L for the platform's attitude measuring system 
are given by substituting the following platform parameters into equation ( 104) : 

R = 2.54 cm 

0 =10.5 deg 

a = 6.0 deg , 

therefore 

L = 2.476 W - 0.0255 . (106) 

Equation (106) was satisfied by letting W equal 25.4 cm and L equal 6.29 cm. 
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Because of the configuration of the attitude measuring systems, certain 
errors are induced in these measurements. These errors are the result of 
measuring angular displacements with a linear measuring device. The nonlinear 
equations that relate the angular displacements 9 of the vehicle with the linear 

displacement (y) of the point of light on the R.T.T. were derived using the 
similar triangle relationship 


y y 

1 L - AR + (W - AR) cos 2 9 


(107) 


Substituting equations (101) and (102) into equation (107) yields 

(W + R)eos0 -R sin 26 

- _ v v 

X (W + R) cos 2 9 + R + L cos 9 - R( 1 + cos 2 9 ) 

V V V 


(108) 


The attitude error that occurs as a result of the transcendental non- 
linearities of equation ( 108) was calculated by substituting the parameters 
associated with the vehicle's attitude measuring system into equation ( 108) . The 
curve "y versus 9^ obtained from equation ( 108) was then linearized to obtain 

a gain (volts/rad) that could be used for the R.T.T. A comparison of the cal- 
culated ■y with respect to the linearized y' for the same 9 indicated a 

maximum error of 2.99 min occurring at 9^= 6. 0 degrees. 

Experimental Maneuvers 


Before any experimental maneuvers were performed with this simulator, 
each of the parameters previously described was measured. Their measured 
values are 

I V = 271 kg-m 2 

I P = 0.257 kg-m 2 

1 S = 9.48 x 10 -4 kg-m 2 

H = 0.627 N-m-s 

K = 130 N-m/rad 


X2 
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c = 


65 


N-m 
rad/ sec 


D = 


0.01395 


N-m 
rad/ sec 


The experiments were performed in the order of their increasing com- 
plexities. The same commands were given to the experimental systems as 
were given to the analog systems studied. Before each experimental maneuver 
was performed, a pretest check of the simulator was made to establish that no 
measureable disturbance torques or unbalanced forces were acting on the system. 

The first experiment was performed using the rigid platform vehicle con- 
figuration and the linear control law. The system's performance curves for a 
step command of 0. 0174 rad are shown by the solid lines in Figure 43. An eval- 
uation of the analog computer and experimental data for all the control system 
configurations will be given later. 

The next series of experiments were performed using the compensated 
platform vehicle configuration and the linear control law. The simulator was 
converted to this configuration by pressurizing the platform's air bearing. The 
control electronics remain unchanged. 

The experimental performance curves for the 0.0174 rad step command 
are shown by the dotted lines in Figure 43. These curves were superimposed 
over the rigid platform system's curves by tracing them from the strip chart 
recordings. The recordings of the maximum attitude maneuver that this system 
can perform (0. 034 rad) are shown in Figure 44. 

Traces of the recorded power and energy curves for the one-degree ma- 
neuver of the rigid and compensated platform systems are shown in Figure 45. 

The complexity of the control system was then increased by adding the 
electronic components that form the nonlinear control law. The simulator was 
converted back into the rigid platform vehicle configuration. The performance 
curves for a step command of 10 degrees are shown by the solid lines in Figure 
46. 


The fourth and most complex control system 
vehicle configuration with the nonlinear control law. 


is the compensated platform 
Again the compensated 
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platform configuration was obtained by pressurizing the platform air bearing. The 
control electronics remain unchanged from the previous experiment. The per- 
formance curves for a step command of 10 degrees are shown as the dotted lines 
superimposed in Figure 46. These performance curves were traced from the 
original recordings. 
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Traces of the recorded power and energy curves for the rigid and com- 
pensated platform systems employing the nonlinear control law are shown in 
Figure 47. 

The last two experiments pertain to the momentum dumping operation. 

The first test was performed using the rigid platform vehicle configuration. The 
control system was commanded to hold an attitude of zero degree while an external 
disturbance was applied to the vehicle. The disturbance torque was applied by 
air currents induced by a variable speed axial flow fan. This disturbance torque 
was adjusted to 0.00678 N-m. The a c torque motor used to simulate the reaction 
control jets was set to develop 14. 15 N for 2.5 sec. The performance curves 
for this operation are shown in Figure 48. 

The simulator was then converted to the compensated platform vehicle 
configuration. The disturbance torque and ac torque motor settlings remained 
unchanged. The performance curve for this platform configuration is shown in 
Figure 49. These experimental curves correspond to Figures 33 and 34. 

DISCUSSION 


A discussion of the passive compensation network and the nonlinear con- 
trol law is presented to draw conclusions regarding the relative merit of these 
systems . 

The passive compensation network is evaluated by correlating the analog 
data of the compensated platform systems with the corresponding data of the 
rigid platform systems. The experimental data are then used to verify the 
differences that were found. 

First, the compensated and rigid platform system employing the linear 
control law is discussed. The one -degree (0. 0174 rad) maneuvers shown in 
Figures 17 and 24 are used as the basis of this discussion. 

The transient characteristics of the vehicle's attitude, as seen by the 
vehicle's position trace in Figure 24, indicate that neither system possesses an 
overshoot. However, the rise time and settling time for the two systems are 
different. The rise time is defined as the time it takes for the vehicle to rotate 
from 10 percent to 90 percent of its commanded maneuver. The measured rise 
time of tiie rigid platform was 8.75 sec, while the compensated platform system 
was 7.4 sec. The settling time was measured to be 21 sec for the rigid platform 
system and 25.80 sec for the compensated platform system. 
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Figure 47. Torque motor's power requirements with nonlinear control law. 
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These transient characteristics agree with the results of the linearized 
analytical study previously performed, where it was found that the time constant 

— — — of the predominant complexed closed-loop poles was greater for the 
^ n 

rigid platform system, while its damping ratio was less. 

The most noticeable difference between the two systems is observed in 
the shape and magnitude of their gimbal velocity traces. The peak gimbal veloc- 
ity of the rigid platform system was measured to be 42 percent greater than the 
compensated platform's gimbal shaft velocity. These traces can also be 
interpreted as being a measure of the gyro-developed control torque, since 
T = 2H<5 cos 6 and 2H in this case is equal to one, so for small gimbal angles 
T . 
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The difference between the gimbal shaft velocities is the result of an 
inherent negative feedback torque that acts on the gimbal shafts. The magnitude 
of this torque is a function of the platform’s or vehicle's angular rate, depending 
on the type of system. Since the platform's mass moment of inertia is 1000 
times less than that of the vehicle's, it develops a much greater initial velocity; 
hence, the compensated system develops the larger feedback torque. 

The peak power required to drive the gimbal shaft torque motor was the 
same for both cases. However, the total consumed energy was greater for the 
compensated platform system, as was shown in Figures 18 and 19. 

The experimental performance curves for these systems are shown in 
Figure 43. The vehicle's position traces indicate an overshoot of 5 percent for 
the compensated platform and 10 percent for the rigid platform system. The 
measured rise time was 10. 25 sec for the compensated platform system and 
11.1 sec for the rigid platform system. The settling times for the rigid and 
compensated platform systems are 35 and 37 sec, respectively. 

The variations between the analog and experimental data can be 
explained by the fact that the gyros used in the experiment had 7 percent less 
angular momentum than those used during the analog study. This infers a 7- 
percent reduction in control torque. 

The recorded power and energy curves for these systems are shown in 
Figure 45. These results correspond to the 0.04-rad maneuver curves of the 
analog study shown in Figure 17 (Sheet 2) . However, the maximum attitude 
maneuver that the experimental systems could perform without exceeding the 
±60 degrees gimbal angle constraint was 0.034 rad. This reduction is also the 
result of the experimental gyros having less angular momentum than those used 
in the analytical study. 

The two traces A0 and A0 were not recorded during any of the com- 
pensated platform system maneuvers. The attitude differential would have to be 
derived from the R.T.T. 's electrical output, and since this device is not per- 
fectly linear, the differential could not be construed as a valid measurement of 
the position and velocity differential. 

The transient response curves obtained from the analog studies for the 
systems employing the nonlinear control law are shown in Figures 29 and 32. 

The gimbal angle and vehicle position traces clearly indicate that the sys- 
tem employing the linear control law is limited by a ±60-degree gimbal angle 
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constraint to a maximum attitude maneuver of ±0 . 04 rad, while the system 
employing the nonlinear control law is capable of performing attitude maneuvers 
as great as 360 degrees. Figures 29 and 32 indicate the 10-degree maneuvers. 

The initial gimbal angle, gimbal velocity, vehicle velocity, and vehicle 
position traces are identical for like systems performing maneuvers of greater 
than one degree. The variable that determines the magnitude of the maneuver 
is the length of time the vehicle is permitted to rotate at its developed constant 
angular velocity. 

The overshoot and settling time of the vehicle's position trace are also 
identical for the same system configurations when performing maneuvers of 
greater than one degree. The overshoot is caused by the nonlinear control law; 
that is, the control law is programed so that active control is not applied to the 
gimbal shaft torque motor until A0 changes from ± to =f . 

The spike that appears on the gimbal velocity traces at the instant of 
command indicates that the developed control torque is an impulse rather than 
a continuous variable torque, as in the case of the linear control law systems; 
however, when the vehicle reaches its commanded attitude, the gimbal velocity 
and control torque become continuous. 

The power and energy required by the torque motor are the same as their 
corresponding linear control law systems for maneuvers of one degree and less. 
For maneuvers of greater than one degree, the power and energy are a constant, 
depending on the type of system. 

The principal differences between the rigid platform and compensated 
platform systems are in the magnitude and duration of the control torque's 
impulse as seen in Figure 32. The magnitude is slightly greater (6 percent) , 
and its duration (0. i sec) is smaller for the rigid platform system. This is 
again because of the larger inherent negative feedback torque created by the 
angular rate of the platform. 

The calculated power and energy curves shown in Figures 30 and 31 
indicate that the peak power required by the torque motor is the same for both 
systems, but the total consumed energy is greater for the compensated platform 
system. This is caused by the voltage (V ) being applied to the torque motor 
for a longer period of time. 



The experimental data for the systems employing the nonlinear control 
law are shown in Figure 46. These results confirm the basic control law 
philosophy and demonstrate the feasibility of physically implementing it into an 
actual flight system. 

The recorded traces of gimbal shaft velocity and vehicle velocity show 
a slower rise time than did the analog data. This is primarily caused by the 
filters that were used in conjunction with the differentiator circuits to derive 
these velocities. 

Another discrepancy is in the magnitude of the developed vehicle velocity. 
The experimental vehicle velocity is less than that predicted by the analog studies 
This condition is created because the gyros are not able to develop as much 
control torque because they have less angular momentum. Because the vehicle's 
velocity is smaller, it required a longer time to acquire its command attitude. 

The traces of the recorded power and energy curves, Figure 47, com- 
pare very well with the analytical values . 

The last area to be discussed concerns the momentum dumping (recaging) 
process . Since this process is independent of the type of control law being used, 
only the system configuration needs to be considered. 

The analog data for the two systems' configurations, shown in Figures 
33 and 34, depict the recaging process for an external disturbance torque of 
0.00678 N-m. 

The dynamic behavior of the system during the period of momentum 
dumping is indicated by the spiked portions of traces. The gimbal velocity is 
again seen to be greater ( 14 . 5 percent) for the rigid platform system . The 
maximum attitude error that results from firing of the thrusters is 0. 00875 
rad for the compensated platform system and 0.008 rad for the rigid platform 
system. The duration of the attitude error for the compensated and rigid plat- 
form systems is 20 and 27 sec, respectively. These differences in system 
characteristics are also obtained from the linearized system transfer functions 
previously presented. 

The recorded experimental data of the recaging process for the rigid and 
compensated platform systems are shown in Figures 48 and 49 . The magnitude 
of the disturbance torque was determined by dividing the maximum angular mo- 
mentum which the gyros were capable of absorbing ( 1. 086 N-m-s) by the time 
it took for the gyro to be precessed to ±60 degrees. Although the experiments 
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were performed consecutively without changing the setting of the disturbance 
torque source, the calculated magnitude of these torques was different (0. 00662 N 
and 0.00715 N-m) . Probably the difference was caused by a change in the 
position of the shielding polyethylene curtains, which were used to keep the air 
currents from the room's ventilation system off the vehicle simulator. 

The experimental recordings confirm the ability of this piece of hard- 
ware to perform the momentum dumping procedure. 

The vehicle's attitude error prior to recaging was greater in both cases. 
This was caused by two things, the lower angular momentum of the experimental 
gyros and a high friction level occurring between 0.4 and 0.5 rad of the gimbal 
shaft. This friction is evidenced by the jagged gimbal angle trace. The magnitude 
and time duration of the vehicle attitude error caused by the firing of the thrusters 
are also greater than those predicted by the analog studies; this again is pri- 
marily a result of the gyro's lower angular momentum. 

The final subject to be discussed concerns the calculated and recorded 
power and energy curves shown previously. These curves were developed to 
ascertain the relative power and energy requirements for the different systems. 
They are not to be mistaken for the minimal power and energy requirements of 
these systems. 

For example, one means of decreasing the energy consumption is to 
replace the gimbal shaft's viscous damping by a negative gimbal shaft rate feed- 
back signal. The gain of the feedback signal is equal to the damping coefficient 
of the eddy current damper. The feedback signal is added to the vehicle's atti- 
tude command. This modification does not affect the dynamic performance of the 
system. 

This method of decreasing the energy requirements was verified by exper- 
imental simulations. The gimbal shaft's rate was derived by differentiating the 
output of the gimbal shaft potentiometer. The eddy current damper was not 
energized during these experiments. 

The performance curves for a step command of 0. 0175 rad applied to the 
modified rigid platform linear control law system are shown in Figure 50. Also, 
the recorded power and energy consumed by the gimbal shaft torque motor during 
this maneuver are shown in Figure 50. 
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When one compares the energy consumption indicated in Figure 50 with 
that of the corresponding unmodified system, Figure 45, the result is a decrease 
of approximately 50 percent in the energy requirements of the modified system. 

This modification can also be applied to the systems employing the 
nonlinear control law. For example, the performance curves for a step command 
of 5 degrees applied to the modified rigid platform nonlinear control law system 
are shown in Figure 51; also, the power and energy curves for this maneuver 
are shown in this figure. 

A comparison of the corresponding unmodified system's energy curve in 
Figure 47 also indicates an approximately 50 percent decrease in the energy 
requirement of the modified systems. 

CONCLUSIONS 


The effects of the passive compensation network on the twin CMG control 
system were obtained from analog studies, experimental simulation, and math- 
ematical solutions to the linearized system transfer functions. The system 
parameters improved by the addition of this compensation network are as follows: 

1. Rise time is faster. 

2. Phase margin is greater; gain margin is equal. 

3. Frequency bandwidth is greater. 

4. Time period of attitude error during recaging is shorter. 

5. Control torques are smaller; gimbal shaft bearing life is greater. 

The system parameters degraded by the addition of this compensation 
network are as follows: 

1. Settling time is increased. 

2. More energy per degree of maneuver is required. 
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3. A larger attitude error during recaging is developed. 

4. Design complexity is increased. 

It should be noted that the margin of improvement of each of the advan- 
tageous characteristics is very small. Therefore, it is concluded that the 
insertion of the passive compensation network into a twin CMG control system is 
not worthy of the additional complexity and effort required to implement such a 
system. 

The advantageous effects of the nonlinear control law on the twin CMG 
system, as determined by the research performed, are as follows. 

1. The magnitude of the attitude maneuvers is not restricted by the 
gimbal angle constraints. 

2. Its attitude acquisition time is faster than a corresponding system 
( same rotor angular momentum) employing the linear control law. 

3. There is no degradation in the attitude (pointing) accuracy. 

4. The control law is easily applicable to the vehicle's mission require- 
ment; that is, the magnitude of the control torque and the linear region of con- 
trol. 


Two features degraded as a result of the nonlinear control law are as 

follows: 


1. The gimbal shaft torque motor has a constant energy requirement 

for any maneuver of greater than ±A0; hence, the energy per degree of maneu- 
ver is excessively large for maneuvers slightly greater than A 0. 

2. The control system's electronics are more complex. 

This control law has been shown to incorporate the best characteristics 
of both the reaction control jet and control moment gyro systems. It is there- 
fore concluded that this nonlinear control law is worthy of consideration for 
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future space vehicles. It is very promising for manned space vehicles used in 
missions requiring fast, accurate, large angle maneuvers. 

George C . Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812 
January 29, 1970 

962-21-01-0000 
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